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SIMU-LC : A LIGHT-CURVE SIMULATOR FOR COROT 
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Abstract 

In order to prepare the analysis of COROT data, it 
has been decided to build a simple tool to simulate the 
expected light curves. This simulation tools takes into ac- 
count both instrumental constraints and astrophysical in- 
puts for the COROT targets. For example, granulation 
and magnetic activity signatures are simulated, as well as 
p modes, with the expected photon noise. However, the 
simulations rely sometimes on simple approach of these 
phenomenons, as the main goal of this tool is to prepare 
the analysis in the case of COROT data and not to per- 
form the most realistic simulations of the different phe- 
nomenons. 

Key words: Data analysis, simulation, p modes, granula- 
tion, magnetic activity, photon noise 



1. Introduction 

Simulating the data that a space instrument like COROT 
will provide might look presomptuous. Indeed, it is cer- 
tainly, when comparing to previous comparable instru- 
ments like IPHIR or GOLF. These two examples show 
that the nominal behaviour of the instrument is not al- 
ways reached, but this does not prevent this instrument 
to provide very interesting data. However, despite some 
technical problems, IPHIR and GOLF yielded a wealth of 
scientific results. Thus, what is the interest of simulating 
COROT data? How close to reality these simualtions will 
get? This might not be the most important fact as the 
preparation of these simulations will help us to prepare 
the analysis of real data and to be ready in case of unex- 
pected technical behaviour of the instrument perturbating 
the data, or unexpected physical behaviour of the targets 
of the instrument. A consequence of that is that the sim- 
ulation tool must include technical and physical aspects, 
making the task even more difficult. These aspects cover: 
photon noise, p modes excitation, granulation signal, stel- 
lar activity signal, orbital perturbations, stellar rotation... 

The software presented here is freely available at: 
www. lesia. obspm. f r/^^corotswg/simulightcurve .html 



2. Photon noise 

The photon noise is certainly the easiest component of the 
noise to simulate... as far as it has the expected behaviour: 
a true white noise with a variance depending on the pho- 
ton counts. The COROT specifications impose a level of 
photon noise of Bq — 0.6 ppm in the amplitude spectrum 
of a star of magitude toq = 5.7 for an observing duration 
of 5 days. The stellar flux for a given magnitude m being 
related to the stellar magnitude by: 

F^Fo io(™o-'")/2-5 (1) 

and knowing that the level of noise varies as the square 
root of the flux, the level of noise is related to the star 
magnitude by: 

S = So 10('"«-")/^ (2) 

As indicated above, this simple relation can model the 
photon noise but not some other photon counting per- 
turbations as slow (periods of hours or more) trends in 
photon countings which might contributes to the low fre- 
quency spectrum of the noise. 

3. Solar-like oscillations 
The solar-like oscillations are stochastically excited and 



simulated here following the recipe of Anderson et al.(1990) 



recalled below. Each solar-like oscillation is a superpo- 
sition of a large number of excited and damped proper 
modes. Each solar-like oscillation can then be decomposed 
as: 

Aj exp[—2iT: Vfjt] exp[—ri{t — tj)]H{t — tj) + c.c (3) 

j 

where Aj is the amplitude at which the mode j with 
proper frequency is excited by turbulent convection, 
tj at which it is excited, ry is the (linear) mode damping 
rate, H is the Heaviside function and "c.c." means com- 
plex conjugate. The Fourier transform of Eq. [3] yields the 
spectrum: 



(4) 

and Aj is a complex number 



1 + 2i{iy - i/o)/r 

whith T = ri/n,U = J2j 
proportional to Aj exp[itj]. 

As the excitations are random, 
Aj has a random phase. As the excitation are very nu- 
merous, according to the central limit princip the real and 



is random and hence 
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imaginary parts of the complex number U are distributed 
according to a normal statitics. Hence the spectrum of the 
oscillations, /(i'), can be simulated by generating the real 
and imaginary parts of the complex number U according 
to a normal statitics. An inverse Fourier transform is next 
applied in order to simulate the oscillations in the time do- 
main. This is the principe of the Anderson et al.(1990) s 
recipe. 

The question is now that is the constraints on the mean 
and variance of U. From Eq.|4]we deduce the power spec- 
trum P(i^) of the oscillation : 

||[/||2 

The mean mode profile < P > is obtained by averaging 
Eq.Oover a larger number of realizations: 

< >^ i + (2(.^.„)/r)» 

where H =< ||?7|p >, which is by definition the variance 
of the complex number U since < U >= 0. We see from 
Eq.[6]that the mean profile has a Lorentzian shape as it is 
- in first approximation - observed. 

According to the Parseval-Plancherel relation, H is 
related to the mean square of the mode intrinsic ampli- 
tude in terms of luminosity fluctuations, Aj^, as (see eg. 



Baudin et al.(2005)) 



H 



where Cg is a visibilty factors, which depends on the degree 
i. Once the mode intrinsic amplitude A/,, the mode line 
width F and Ci are given, we have then a constraint on 
the variance of \J . 

Note that the derivation of Eqs.[3][7] assumes that the 
the mode life-time (^ ?7^^) is much shorter than the dura- 
tion of the observation (Tq), that is Tq 77 ^ 1. For life-times 
much longer than Tq , different expressions are derived but 
the principe of the simulation remains the same. 

For sake of simplicity, we use the adiabatic assumption 
formulated by Kjeldsen fc Bedding(1995), to deduce the 
maximum of from the maximum of the root mean 
square mode velocity Vmax according to: 



(^l)© max \r 



?0,cff 



(8) 



0,max V Pcff 

where Toff is the effective temperature and the symbol © 
refers to quantities related to the Sun. We take for the 
Sun the rms values (^L)0,max — 4 ppm (see table 2 in 
Kjeldsen & Bedding(1995) and more recently Barban et al. 
and Vo^max — 27 cm/s according to Chaplin et al.(1998) s 
seismic observations. 

In turn the root mean square of the mode intrinsic 
velocity, V, is related to the rate V at which energy is 
injected into the mode by turbulent convection and the 
mode damping rate ry = tt F as: 



where Ai = I /^"^{h) is the mode mass, / the mode inertia, 
the mode radial displacement and h the height above 
the photosphere where oscillations are measured (we con- 
sider /i = 0) . The way the quantities involved in Eq.[9] are 
model is explained in the next two Sects. 

3.1. Mode excitation and damping rates 
Theoretical mode damping rates are obtained from the ta- 



bles calculated by Houdek et al.(1999) These calculations 
rely on the non-local and time-dependent formulation of 
the convec tion of |Gough(1977)j |Gough(1976")| (for more 



(5) details, see [Houdek et al.(1999)P 



The computation of the excitation rates V is performed 



according to the model of stochast ic excitation ofjSamadi fc Goupil 
The calculations assume 



Samadi et al. (2003a) 



as m 

a Lorenzian function for modelling the convective eddy 
time-correlations. Furthermore, the characteritic wavenum- 
ber fco involved in the theory is assumed to be constant ac- 



cording to the simplification proposed by Samadi et al. (2003b) 
Its value is related to the value /cq,© ~ 3.6 Mm~ inferred 



m 

as: fc( 



Samadi et al. (2003b) from a 3D simulation of the Sun 



;o = ko^Q KqI K where A is the mixing-length evalu- 
ated at the layer the convective velocity is maximum. 



3.2. Mode mass and mode frequency calculation 



(7) The solar model we consider is calculated with the CE- 



SAM code (Morel(1997) ) and appropriate input physics, 



described in details in Lebreton et al.(1999) In particular. 



convection is modelled according to the classical mixing- 
length theory Bohm - Vitense(1958)[ iereafter MLT with 
a mixing-length A = aHp where Hp is the pressure scale 
height and ac is the mixing-length parameter. In contrast 



with Lebreton et al.(1999), the atmosphere is restored from 
the Eddington classical gray atmosphere and microscopic 
diffusion is included according to the simplified formalism 



of ,Michaud fc Proffitt(1993)[ The calibration of the solar 
model, in luminosity and radius for an age of 4.65 Gyr, 
fixes the initial helium content Y — 0.2751, metallicity 
Z = 0.0196 and the MLT parameter = 1.76. 

The oscillation eigenfunctions and hence the mode ma- 
sses, M, in Eq.iniare calculated with the adiabatic ADIPLS 



pulsation code ( Christensen-Dalsgaard et al.(1996)). 



1/2 = 



4. Stellar granulation signal 

tion can be considered as noise if the aim of the ob- 
servation is the stellar oscillations, but it carries some in- 
formation about the physics of the star, and very valuable 
information as the convection is a very badly described 
phenomenon in stellar modelling. Thus, we call it a signal 
and not a noise. 

However, granulation can be described only with a sta- 
tistical approach. Moreover, its contribution in the Fourier 
domain is not independent of frequency: it will contribute 
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imulation/AII components 
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Figure 1. Left: Fourier spectrum of the output of the light curve simulator described here, in the case of the Sun; Right: 
Fourier spectra of observation of the Sun obtained with the photometer LOI onboard SoHO (the high peaks at low 
frequency are daily aliases due to SoHO synoptic observations performed every 24 hours). This comparison shows a 
reasonnably good agreement, except maybe for the activity signature 



more at low frequency. A common description is to con- 
sider that signal S" as a random signal with some memory: 
then, its autocorrelation function is: 



ACFs{t) = 



\t\/r 



(10) 



where A is the amplitude of this signal (to be related to its 
variance) and r a characteristic time. Knowing that the 
Fourier transform of the autocorrelation function is the 
squared modulus of the function S, one has: 



15(^)11 



2A^ 



(11) 



1 + (27rTi/)2 ■ 

Thus, the Fourier spectrum of the granulation is modelled 
as a Lorentzian which, if it is summed over all frequencies, 
yields the variance of the signal and then its relation with 
A: 



= J \\S{iy)\\ ^ (7 = A/V2. 



(12) 



So, knowing the intensity standard deviation due to gran- 
ulation and its characteristic time, it is possible to model 



its Fourier spectrum (and considering its phase is random, 
it is possible to model the corresponding time series by an 
inverse Fourier transform) . This is in fact the approach de- 
velopped by Harvey(1985) to model the solar noise spec- 



trum, including granulation, mesogranulation and super- 
granulation. It can be used to describe any noise with 
memory (some electronic noise for example). 

In the case of granulation, the required parameters are 
estimated from theory of convection. The eddy size Acddy 
is assumed to be equal to A = a Hp where a is the mixing- 
length parameter and Hp the pressure scale height. The 
eddy overturn time is related to the eddy convective ve- 
locity, V, by Teddy = Aeddy/w and is considered as the char- 
acteristic time of granulation r. The number of number of 
eddies seen at the star surface, A^eddy, is 2 7r(i?,/Aeddy)^ 
where i?* is the star radius. This relation ignores of course 
that the medium is higly anisotropic. According to the 
mixing-length theory (see eg. Cox(1968)|, the eddies con- 
trast (border/center of the granule), (SL/L) eddy, can be 
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related to the difference between the temperature gradi- 
ent of the eddy, V and that of the surrounding medium, 
V, according to the relation: 

(<5L/L)eddy = (13) 

In turn, Eq. [T3]can be reduced to: 

(<5i/L)eddy-(4/9)-A^i (14) 

(1 - A) 7 

where 7 is the convection efficiency and A the ratio be- 
tween the convective flux and the total flux of energy. 
This relation is finally calibrated in order to match the 
solar constraints, and the intensity standard deviation for 
the whole observed disc is 

a = (<5L/L)eddy/\/A^cddy (15) 

All the quantities are obtained on ID stellar models com- 
puted as explained in Sect. [31 

5. Stellar activity signal 

The stellar magnetic activity will induce intensity varia- 
tion in time, mainly due to the presence of starspots cross- 
ing the observed disc because of the rotation. However, 
some other sources of intensity variations are expected 
(flares for example). 

A first approach is to consider the intensity variations 
as described by the same way than granulation. Knowing 
the standard deviation and the characteristic time of in- 
tensity variations allows to build the Fourier spectrum of 
these variations as a Lorentzian. This approach has been 



used for example by Aigrain et al. (2004) The difficulty in 



the case of magnetic activity is that there is no theoreti- 
cal description of the phenomenon. Thus, the parameters 
describing it are empirically derived. 
The characteristic time is taken as the period of rotation 
of the star, or, if the rotation is slow, the instrinsic life- 
time of a spot. The latter is arbitrarily chosen as the one 
in the solar case, computed by Aigrain et al.(2004) The 



rotation period is computed from an empirical law involv- 
ing the age and the B — V color index of the star described 
in the same reference. 

The standard deviation of intensity variations due to mag- 
netic activity is also derived from an empirical law. How- 
ever, this law involves the Rossby number Rq (ratio of ro- 
tation period Prot and the overturn convective time at the 
base of the convection zone Tbcz ) ■ This number can be de- 
rived empirically Aigrain et al.(2004)[ but in the present 
simulations, we derive Tbcz from models. Then we use em- 
pirical laws to estimate a. 

Another approach is to simulate in the intensity time 
series the influence of individual spots, to estimate their 



number and contrast (as for example Lanza et al.(2004)) 



The expected result in the Fourier spectrum should be 
similar to the first approach, but this detailed simulation 
should allow for rotation measurements. This approach 
will be included in a further version of our simulation soft- 
ware. 
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Figure 2. Fourier spectra of the output of the simulator 
described here, in the case of a Sun seen with a magnitude 
m — 6, showing that the p modes are easily detected, as 
well as granulation and activity 



6. Examples 

A first example of the output of this software is shown in 
FiglH a simulation of a Sun with a magnitude m = (in 
order to have a very weak phoon noise) is compared with a 
spectrum from the LOT photometer ( Appourchaux et al.(1997) I 
onboard SoHO. The agreement between simulation and 
ovservation is good, except maybe at very low frequency: 
the activity component of the signal is overestimated in 
the simulation (this is explained as the Sun does not fit 
very well the empirical law used for activity estimation). 
The following example is a sun with a magnitude m = 6, 
showing that both p modes and activity should be visible 
with COROT for such a star (FigH]). Another example is 
shown in Figl3]for a star of 1.5 solar mass and an age of 
2.4 Gyr: again, p modes and activity are detectable. 



7. Conclusion 

This simulator software will continue to evolve with time. 
As indicated above, intensity modulation due to starspots 
will be included, as well as other stellar or instrumental 
signals, as for example instrumental perturbations due to 
orbital vraiations. Moreover, this effort of simulation will 
not end with the delivery of first data but will be contin- 
ued after that. The comparison with real data will allow 
to check for the validity of physical hypothesis used to 
simulate the different signals of astrophysics origin in the 
data. This shoud bring a great amount of information on 
our knowmedge of these often not well known phenomena, 
which stellar simulation is often derived from the solar 
case. In parallel, the simulation of instrumental compo- 
nents of the signal will be improved to help the interpre- 
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Figure 3. Same than Fig. d for a star with a mass M = 
1.5-Mq and a magnitude m = 8, showing also p modes, 
granulation and activity signatures 
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